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ABSTRACT 


The accuracy with which the earth's gravitational field can he ■ 
defined from an analysis of satellite data is strongly dependent 
upon the satellite constellation used. The constellation selected 
must satisfy several criteria. The objective of this study was to' 
produce quantitative information on the effects of the geopotential 
harmonics on the motion of satellites for a large variety of orbits. 
This data can then.be used to investigate satellite constellations 
which will yield useful geodetic information. 


-iii- 



THE JOHNS HOPKINS UNIVERSITY 

APPLIED PHYSICS LABORATORY 

SILVER SPRING MARYLAND 


■eniMG PAGE BLANK NOT FILMED. 


TABLE OF CONTENTS 


Page 

I. INTRODUCTION 1 

II. DEFINITIONS AND NOTATIONS -3 

III. DISCUSSION OF RESULTS 5 

Acknowledgments . . . . . . 13 

References ........ 15 

APPENDIX A - Expansion of the Gravitational Potential 17 

APPENDIX B - Derivation of H, C r L Equations . . 25 

APPENDIX C - Graphs of Inclination Function . 35 

APPENDIX D - Graphs of Inclination Function F^ . yy 

APPENDIX E - Graphs of Eccentricity Function . . 119 


-v- 



THE JOHNS HOPKWS UNiVERSiTY 

APPLIED PHYSICS LABORATORY 

SlUVER SPRING M*«Y1>H0 


I. UTOQDUCTION 

In future satellite geodesy programs to improve on the current models 
of the geopotential, additional satellite orbits will be needed. Some of 
the present models of the earth’s gravity potential contain over 200 
coefficients. It is fairly obvious that as more terms are added to the 
model the problem of separating out their individual effects becomes 
increasingly difficult. Separation of these harmonic effects is strongly 
dependent on the orbital characteristics of the satellites used in the 
analysis. Therefore serious consideration should be given to the problem 
of defining a satellite constellation to be utilized in a satellite geodesy 
program. 

There are several factors to be considered in selecting this satellite 
constellation. Since the major expense of a satellite geodesy program is in 
the cost of the satellites and their launch vehicles, it is important that 
the constellation consist of the minimum number of satellites that is 
needed to adequately determine the geopotential model. This implies that 
maximum use of existing satellite data should be made. Additional satellites 
therefore should be selected to complement the existing data, which are of 
sufficient quantity and quality. 

The main objective of this paper is to provide quantitative information 
on the contributions of the harmonics in the earth's gravitational field to 
the motion of a near-earth satellite. The harmonics through 20th degree 
and order are considered. This provides information necessary in selecting 
a constellation to be used in a satellite geodesy program. Most of the 
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analyses are contained in the Appendices. The notation used is given in 
section II and the main results are discussed in section 'HI. 
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II. DEFINITIONS AND NOTATIONS 

x,y,z = inertial cartesian components of the satellite 

r,cp,\ = spherical radius, latitude, and longitude bf the satellite 

H,C,L = components in a -right handed coordinate system with axes in the 
three directions 

1 - radius vector from Earth* s C.G. to instantaneous position 

of the satellite 

2 - negative angular momentum .vector direction 

3 - component of satellite velocity perpendicular to radius 

vector and lying in the orbital, plane of the theoretically 
computed orbit 

a = semi-major axis of the satellite orbit 

e = eccentricity 

i - inclination 

to = argument of perigee 

Cl = longitude of ascending node 

f = true anomaly 

p = f + to 

M = mean anomaly 

E = eccentric anomaly 

n = mean motion 

Aq = right ascension of Greenwich 

U(r,cp,\) = potential function 

p - gravitational constant for the Earth (= 398601.5 fcnr/sec j 

P^ m (sincp) = associated Legendre polynomial of degree i and order m 

= unnormalized coefficients associated with the spherical harmonics 
of degree and order ( £,m) • 

N^ = normalization used for spherical harmonics 
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N 




U-m)> (2JJ+1) (2-5 0lil )l l/2 

(j£+m ) \ 


om 


fo for in / 0 
*1 for m = 0 


R 


near equatorial radius of the earth ( = 6378 -. 144 km) 


Zero subscripts are used to designate parameters associated with a reference 
orbit . 
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III. DISCUSSION OF RESULTS 

In, selecting a satellite constellation ’to be used in a program to deter- 
mine a model of the earth's gravity field, it is important to understand 
the physical effects -each harmonic in the geopotential expansion will have 
on the satellite's motion. Studies have already -been made on using resonant 
orbits for a geodetic program [.1,2], Resonance occurs when the satellite 
orbital frequency becomes commensurate with an integer multiple of the 
earth’s sidereal rate. This results in greatly amplified motions of the 
satellite due to those harmonics which are in resonance, thus allowing for 
a more accurate determination of their coefficients. However this approach 
would require a considerably larger satellite constellation,- since a given 
orbit can be selected to resonate with only a relatively small subset of 
the geopotential harmonics. As the major cost in a satellite geodesy 
program is in the satellites and their launch vehicles, the use of resonant 
orbits becomes undesirable, financially. 

It is important, however, in future attempts to improve the geo- 
potential model to incorporate data from satellites which are more sensi- 
tive to the effects of the earth’ s gravity field., particularly to the higher 
order harmonics. This can be achieved through the use of very low altitude 
orbits. Until recently, these orbits were not useful because of the 
unpredictableness of atmospheric drag which contaminated the data. Con- 
siderable effort has been put into the development of a thrust package to 
be incorporated into. the design of a satellite which would counteract the 
effects of all non-gravitational forces acting on the satellite [ 3 ]. This 
device woulcl make feasible the use of low-altitude orbits as well as simplify 
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the analysis required to reduce the tracking data received. Since the effects 
of radiation are also eliminated, "the inclusion of this package on satellites 
at all altitudes would he desirable. 

To study the effects the harmonics in the geopotential expansion will 
have on the motion of a satellite, their contributions are derived in a 
coordinate system moving in inertial space. The origin is located at the 
satellite as defined by a reference orbit. The vector difference between the 
true position and the position as defined by the reference orbit is parti- 
tioned into the H, C, L components ©f this system. (Definitions are given 
in section II.) Contributions of the geopotential harmonics in this system 
are derived in Appendix B. The final results given in equations (B.22) - 
(B.21}-) are restated here for discussion purposes. 


H 


e 


(l-e 2 ) 


jUt 2 


R JL + 2 . — , r— > \i!.j 

(-y) Y V I (i) F (e) 

a Z_. 0 or. 


A 


(0,0) (H) 


p=G k= - 00 


Jimp 




5 A,m,p,k 


(B.22) 
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From the definition of in Appendix A, equation (A.2l), we made the 

substitution 


1 + 


1 - 



.(1 - e 2 ) 


- 1/2 


in the above expressions. The S-functions are defined as 

( X-m)even 


r- l-sh 

- <- x ) 1 2 J ' 


77 m 

C X 


o' m 

“ S X 


c m 
S X 


77 m 


cos [(X-2p+k)M + (X-2p)o + mCd-A^,)] 


(X-m)odd 

(X-m)even 

( X-m)'odd 


sin [(X-2p+k)M + (X-2p)w + m(d-A^) ] 


For a given satellite orbit (to the accuracy of the solution's), the time 

dependence of the solutions are contained in these S-functions. They 

encompass all possible frequencies that can be generated by a harmonic of 

degree X , and order m , The summing index k denotes those frequencies 

k. 

generated by terms of order e and higher. For each frequency the inclina- 
tion and eccentricity dependence can be factored into. the I and F 

¥ 

functions, respectively. The I-funetion is defined in equations (A.l8) 
and (A. 19 ) The F- functions are -given in equations (B.26)-(B.28) . 

Because of their complexity they will not be rewritten here. The frequency 
contributions are imbedded in the F- functions. From their given form, it 
is apparent that the amplitude of the contribution of each oscillation is 
inversely proportional to its corresponding frequency. Therefore, in general, 


(A.28) 
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the low frequency terms associated' with each harmonic will he the -ma^or 
contributors to the satellite’s motion. Thus • ty examining in detail the 
amplitudes of these low frequency terms, significant insight into the effects 
of the geopotential harmonics on the satellite’s motion, can be obtained. 
With the given partitioning of the inclination and eccentricity effects, 
one can investigate a much larger variety of orbits. 

The possible frequencies that arise for a given and m can be 
written as 


(£-2p+k) M + (4-2p) ® + m (Q-A g ) (l) 

As a general rule, the larger amplitudes are associated with the low frequency 
terms which occur when k = 0 . This follows from the fact that k incor- 
porates those frequencies which result from- the eccentricity -effects- of the- 
harmonics, (it should be noted that these low frequencies associated with 
the harmonics of a given older m and of the same parity in degree 1 
will be identical. This is the main reason why it. is difficult to separate 
any two of these -harmonics using data from only a single satellite.) 

In the H and L equations, (B.22) and (B.24) , the low frequencies 
investigated were for values of J&-2 p equa-l to 0 and 1 . The inclina- 
tion dependences in these components are identical and the graphs associated 
with this function are given in. Appendix C where 

(0,9) 

F 2 = I (i) (2) 

Mp 
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Each graph contains one or more plots of F^ versus inclination for a 
specified value of m for the range 0°^ i S l80°. Each curve corresponds 
to a different, value of & and the set of curves in each graph consists of 
those H values which yield the same low frequency terms for the given m . 
This grouping is helpful in selecting orbital inclinations which will give 
good separation of the harmonics represented. The graphs indicate that the 
contributions of the harmonics of order 13 and higher for satellites with 
orbital inclinations less than 40° and greater than 150° will be negligible. 

The inclination dependence in the C component is given by 

( 1 , 0 ) ( 0 , 1 ) 

F^ = 2 cosi I (i) - m sini I (i) (3) 

Xmp Mp 

Similar graphs of this function are presented in Appendix D for values of 
Ji-2p~l equal to 0 and 1 . Since the inclination functions are independent 
of k , the values generated are applicable for the frequencies corresponding 
to all values of k . 

To obtain meaningful information on the effects of varying eccentricity 
it is necessary to consider simultaneously the semi-major axis contributions. 
This is accomplished by requiring that the perigee altitude remain constant 
for all eccentricity values considered. You can then express a as a 
function of e and the fixed perigee altitude. For computational purposes 
we arbitrarily fixed the perigee altitude at 638 km so. that 
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By also defining 


n-A 


G = (A-2p=)“ + m (-—) 
h n 


Q-A, 


G = (A-2p-l) -r- + m (-r-^) 
n n 


w 


• — ' ♦ * 

the explicit dependence of the functions' on m, (*>_, and Ag are contained 
in a. single parameter. The final form of the functions plotted for each 
component is given by 


F(H) - 


£+2 

l-e\ 2 


Ti.i) 


f— ') 
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F 


(H) 

(e) 

-^G^p^k 
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F(C) = 


( 1 . 1 ) 
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X+2 

(i£) 2 

^ T Me* * 


(C) 

F (e) 

4,G 2 ,p,k 


X+2 


(L) 


F(L) =-^-^2 (~) 2 P' '(e) 

(1.1) ^°l^ k 


( 6 ) 


(T) 


In Appendix E these functions are plotted versus .e for the range -0 < e s 0.5* 
The parameters G^ and G^ take on values in the range -0.1 ^ G^, Gg ^ -1.5 
in increments of 0.2. Gg appears only in the C component, G^ is con- 
tained in the H and L expressions. The same values of t-2p -were used 
with k set equal to zero. Each graph contains a set of curves associated 
with a given value of & and each curve is generated for a constant value 
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of (K . (i = 1 or 2, depending on which function is being plotted.) In 
plotting the functions, a distinct symbol was used for each value - of . 

This relationship is presented in Table 1. 

The use of the G. parameters increases the amount of information 
contained in each graph. However , to relate these graphs to a particular 
harmonic requires an additional computation to obtain the appropriate values 
of G^ . Equations ( 4 ) give the formulas needed. The range of values used 
for G^ should cover the majority of orbits of interest for investigating 
useful satellite constellations for geodetic studies. The effects of a 
different value of the semi -major axis can be incorporated by multiplying 
the approriate number from the graphs by 


K 


a 


1.1 R 

e 

a(l-e) 


4+2 


( 8 ) 


The magnitudes of the functions decrease monotonically as the eccentricity 
increases, for a majority of the cases plotted. We also evaluated the functions 
for k = -1 . These results yielded increasing magnitudes as e increased, 
however their maximum values were much smaller. These computations are not 
included in this report. 

The information contained in this report is by no means sufficient to 
base the selection of a satellite constellation to be used in a geodesy 
program. There are many other factors that must be considered, for example, 
as mentioned earlier, satellite and launch vehicle costs, and usefulness of 
existing data influence the selection. Equally important is the coverage 
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Table 1. 

Symbols used to- designate the 
value of Gp and G2 corres 
ponding, to eaph .cinnve in ' 
Appendix E. 

G 1 ' G 2 


Symbol 

-O.l 


0 

-0.3 


X 

-0.5 


it 

-0-T 


+ 

1 

O 

VO 


X 

-i.i 


A 

-1.3 


O 

-1.5 


E3 
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available and quality of data obtainable for a given orbit. However ^ proper 
analysis of these data can yield valuable information for purposes of orbit 
selection. 
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PRECEDING PAGE BLANK NQT FSL^Vp/ 
APPENDIX A 

Expansion of the Gravitational Potential 


The earth 1 s gravitational potential is defined in terms of* a series 
expansion in spherical harmonics. 




cos m(\ - A g )+ s/ sin "m(X 


The notation used is given in section II. In examining the gravitational effects 
on the motion of a near-earth satellite, a normalization factor is introduced 
for the surface spherical harmonics which is reflected in the determined values 
of their corresponding coefficients (C^ m , ). We have chosen the normal- 

ization which makes the integral of the square of each harmonic over the 
surface of a sphere to be Ltt [ L] . This is given by 


where 


(A-m)> (2A+1) (2-6 Qm ) 


0 for m ^ 0 


1 for m = 0 


a m _ 


1 r m 
1 


The mean equatorial radius of the earth, , is used to define the unit of 
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length, and the unit of time is (2tt) P times the Schuler period. This 
time transformation is given by 




(A. 4 ) 


The gravitational constant, p , expressed in these units is unity. Incorporating 
the above changes, equation (A.l) becomes 


E j& 


& a 

U = ~ (•“) {1 + X! P / ( sinc P)^/ 1 cos + s/ 1 sin m (^-A G )l] 

i,=2 m=0 

(A. 5) 


which is partitioned further into a central term and perturbations about this 
term: 

u - f £> ♦ V 
e 


(A. 6) 


Transforming to the kepler elements, V can be expanded into a form more useful 
for a perturbational analysis. For a harmonic of given degree and order, (A,m) , 
we define the complex function 




+ j v 


m 

JL 


where 


u 


m 


Jt> 


v, 


m 


P^ 11 (sin 9) cos m(X-A G ) 
P^ 331 (sin 9) sin m(l-A G ) 


(A. 7 ) 


(A.8) 


such that 
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(A.9) 


From spherical geometry we have that 


sin cp = sin i sin P 

< cos cp cos(X-Q) = cos P (A. 10) 

cos cp sin(X-fi) = cos i sin P 


Also, the Legendre polynomials can be rewritten as 


P^ m (sin cp) = cos m cp T^ m (sin cp) 


(A. 11) 


where 


m m / . - \ 1 

T ji ( sin <P) = £1 


rill] 

L 2 J 

.1 

t=0 


(-1J (2A-2t)! [sin i sin p] 

t l ( l- t ) I ( A-m-2t ) \ 


j£-m-2t 


(A. 12) 


and 


[A=] - / 


A-m 

2 

#~m~l 


for JL-m even 


for £-m odd 


Incorporating equations (A.1G) and (A.ll) into (A. 7 ) yields 


= T™ (sin cp) [cos p + j cos i sin P] ^ 


(A. 13) 
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Wow, . replacing 


f sin ^ - (-gi) [e jP - e'&) 


cos P = (|) [e^ + e' JP ] 


into eq. (A. 13 ) and expanding, we obtain 


U-m] 
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-3 I 


j£-m~ 2 t 

( -1 ( 2 A- 2 t)l sin i J ^)X-m- 2 t 


t! (A-t)-l U-m- 2 t)V 


t =0 


m A-m- 2 t 

I I 

s =0 c =0 


2 (m-s) 2 s ■ , 

(-l)- c cos, i sin i (“)( “1 ) exp. j{[ £- 2 (t+s+c,) ]p 

.2 2 


+ ni(Q-A G )3 


The above equation can be rewritten so as ‘to collect the coefficients of each 
argument of the function exp j {' [ A-2 ( t +s +c ) ] (3tm( fi- ) } . - Let 


p = t+s+c , 


where 0 ^ p 1 £ H . In adding the summation on p we can eliminate one of 
the other summations. In this case we set s = p-t-c and eliminate the 
summation on s . Finally, 


(A. 14) 


(A. 15) 


(A. 16 ) 
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(A. 18) 
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VI 
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(A. 19) 


The additional parameters g and h ' are introduced into the definition in 
order to facilitate the computation of the forces derived from the potential. 

We now need to- expand the functions of r and f as a series in the 
mean anomaly, M . Utilizing the derivations given in [ 5 ] we can obtain 


r 11 e*^ - 
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<1 = A - 2p 
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a(a-l)» » ’(a--k+l) 
kl 


Equation (A.20) is valid for all integer values of q and k . Eq. (A. 24) 
is Just the definition of the Bessel functions of the first kind and eq. (A. 25) 
is one form of the hypergeometric series. With these expressions, we. can now 
write the perturbing potential for a given 4 and m . 
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(A. 26) 
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(A. 28) 
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preceding 

APEEHDIX B 

■Derivation of H, C, L Equations 


The effects of the geopotential harmonics on the motion of a near-earth 
satellite are derived in a coordinate systemt moving in inertial space. The 
origin of this system is located at the .satellite’s position as defined hy 
a reference orbit. We denote the true satellite position vector by r 
and use a zero subscript to characterize ; those orbital parameters defined 
by the reference orbit. We now define 


Ar = r - r = Hh+C£+Lje 


where 

r 


A 

h = 


\ C = 73 ? 


dr 

o — 
df3 0 ° 


dp O 


A 

A 


dr 


r x (~ x r ) 


o dp 


dr 

(^o 2 x rJ 


(B.l) 


(B.2) 


The H, C, „L - space, so- defined, is an orthogonal right-handed coordinate 
system. It is assumed that errors in defining the reference orbit are due 
to errors in the geopotential harmonic coefficients, denoted by AC™ and 
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AS™ , and are of the order 0 (10 - ^) . A first order perturbation theory 
is used, therefore terms which contain C™ AC™ , S™ AC™ , etc. are omitted. 

The radius vector in inertial cartesian coordinates is 

r =<xi+yj+zk ( B *3) 


Setting x = x + Ax , etc. in eq.. (B.3) and substituting into (B.l) gives 


A A A A A A 

Hh+C£ + Lj& = Axi + Ay«j + A z k 


(B.4) 


or 


I A 

c 


\ L i 


/a a 

'h * i 

A- A 

£ • i 

A A 

A ■ i 


A A 

h • j 

A A 

3 

A A 

A • 3 


A A 

h • k 

A A 

k 

A A 

A • k 


/ Ax\ 

.Ay 

\ AZ / 


(B.5) 


In spherical coordinates 


r 

X 
z 


r cos cp sin X 
r cos <p sin X 
r sin cp 


(B.6) 


Setting r = r Q + Ar , etc. into (B.6), one can obtain 


/ cos cp cos X -r sin © cos X„ ~r cos <p sin X. \ A r\ 

I T /Vi (T) A 1 {T\ '(Cl 


cos cp sin X 
Y o o 

sin cp 


O 'O o 

-r sin cp sin X 

ft M r. 


V cos cp • 
o 


r cos cp cos \ 
o^oo 


-0 


/ 


Acp 

\ 4X / 


(B.7) 
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Substituting (B.7) into (B.5) and simplifying, .yields 


- c 

J 

l L j 

\ 


0 

cos i Q 

- r o cos .cp Q o 

sin i cos P 
r o © 

cos 9 


0 ^ 


M 
' ' 

sin i cos P 
0 0 


A9 

1 r cos i I 


W 

O O/ 


or taking the inverse of the matrix 


M 

. Aq> 

W 


h 


\ 


0 

-cos i 


r cos cp' 

sin i 0 

r cos cp 
o T o 


° . \ 


sin i cos p 

O _ Q 

r cos cp 
o o 

cos i 
o 

r ' coi^ qp 

o ^ 


1 


/ h \ 

C: 

w 


(B.8) 


(B.9) 


At this' point in the analysis there are two- ways in which one can 
proceed. The first approach results in a set of differential equations 
relating, the H', C r L coordinates to the geopotential harmonics. 'To obtain 
these equations we start with the equations of motion of a satellite in 
spherical coordinates. 


r 






av 

ax 


(B. 10 ) 
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In the above we again make the substitution 


v = r + Ar 

o 

9 = 9 0 + Aq> 

X = X Q + AX 

V = V + AV 

^ o 


(B.ll) 


a nd expand' to first order in the differential elements. The'zer©^ order term 

in the equations cancel since they are the defining equations of our refer- 

■> 

ence orbit. The resulting first order equations are 


Ar - 2r cp (Acp) - 2r cos 2 cp (xO(AX) - ['.(<p) 2 + cos 2 cp (X _) 2 + Ar 

O t O T O O O 0 O , D 

O 

. 0 /-x2 A SAV 

+ r-em-lStM 4¥--g- 

o 


< (A<p) + 2r Q (qn )Ar} + r 0 2 (X Q ) sin-2cp 0 (AX> + r 0 (X/) 2 «ln 2<p 0 Ar 

2 ,: s 2 o a ^av 
+ r o (X o } C0S S<P o A<P = W 

O 

k. rrr {r 2 cos 2 cp (AX) - r 2 sin 2cp (X ) Acp + 2r (.X' ) cos 2 cp Ar] - -tr— (B.12) 

dt c o Y o o o T o' o T © 3X 0 


In these equations the reference, orbit can be approximated by an ellipse. 
Therefore substituting .equation (B.9) into .('B.12) and simplifying gives 
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H - [(p) 2 + H - 2 r (p) jt (— ) = 

r 3 ' dt v r Q ' dr Q 


L - — U + 2(p ) H + 2 [— P +> M H = ™ |~ 
r v o' r o ) o r SP 

0 O' GO 


c + [(p) 2 


“] c ~- 

r o r o Sinp o 


dAV 

51 o 


(B.13) 


Solutions to these equations are best obtained via numerical integration 
techniques. 

The problem of solving coupled second order differential equations 
however can be avoided. This alternate approach provides solutions better 
suited for the objectives of this paper. We first define the H, C, L 
coordinates as a function of differential kepler elements.- The satellite 
radius as a function of the osculating kepler elements is given by 


r = a (l-e cos S3) (B.llj-) 

where the eccentric anomaly, E , is related to M , the mean anomaly, through 
Kepler’ s equation 

E - e sin E = M (B.15) 


These equations yield 


Ar 


= ~ Aa 


a 


a e sin f 
(l-e 2 ) 



AM - a cos f Ae 
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-From equations (A.10) in Appendix A and 


cos f = 


cos E - e 
1 - e .cos E ; 


(B.17) 


which relates ihe true anomaly* f to the eccentric anomaly } we obtain 


^ Acp = sin i cos_P + cos i sinjg M 


cos cp 


cos cp 


/ AX = 


k 1 cos i A q 
AO + g - Ap 

cos cp 


cos p Sin cp .. 

2 

cos cp 


AP = Au> 


, jx±£±!l 


AM * Siaf [*0^-1 + !] Ae 
2\ r 


(1-erf) 


(B.18) 


Substitution of equations (B.l6) and (B.l8) into (B.8) yields 


'"h - £ ia + “ 8 w/ AM - a cos f Ae 

a [l-e 2 ] 1/2 


\ C = r [sin i cos P AO - sin P Ai] 


L = r [Au> + cos i AO 


, a^l-e 2 ) 1 / 2 


AM] 


(B.19) 


These differential kepler elements also ■ represent the contributions of the 
errors in the force model used .to generate the reference orbit Their 
coefficients in the above equations can be evaluated using a set of mean 
kepler elements. The differential kepler elements are related to errors 
in the geopotential’ model through a perturbation on Lagrange’s planetary 
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equations of motion. To first order they can he • written as 


Aa - 


Ae 


n.(l-e ) 


sa/s 


[e sin f AR + 


i(l-e ) 


AS] 


. (1 -e g ) 


2x1/2 


na 


[sin f AR + (cos E + cos f) AS] 




Ax = 


sin i .AQ - 


Att) = 


V. 


AM 


1 


naCl-e 2 ) 1 / 2 - 8 


- cos AW 


nafl-e 2 ) 1 / 2 - 8 

q-e 2 ) 1 / 2 


sin P AW 


nae 


[-cos f AR + ( + 1) sin f .AS] - cos i AQ 

a(l~e 2 ) 


1 1 Aa - ^ AR - (1-e 2 ) 1 / 2 [ip ,+ cos i AQ] 
na 


where 


r 


AR 


5AV 

dr 


{ .AS 


i MZ 

r SS- 




AW 


i m 

r sin£ di 


Using the results obtained in, Appendix A , equations ("B.20) are integrated 
and then substituted into (B. 19 ). After collecting terms with the 
frequency., the final expressions become 


(B.20) 


(B.21) 


same 



l+£ 2 1+2 R 1+2 A Z 


An " -2 _ 2 - i 

n R e 1_P l 


£z±-) <JL) y y i (o ' o) p <«> (o s 

2 ' a 2_. 2_. Anp AjG^fcV A,m,p,k 

p=0 &,= ■* 


l+P 2 1+2 R 1+2 A “ 


C, = “#—2 ( (~) Y V [2 cos i - m sin i F (e) S. . . 

An -2 R 2 2' (_> Anp Anp £,G£,p,k A-l,m,p,k (B.23) 


e - "1 


prO ]C = -« 


1+P 2 X+2 R £+2 JL A 


l a-r^x.) &> y y iMjm 'Ds 

n^ 2 l-6 2 L L X Anp P A,G[^,k v_1) A>,p,k 


e 


p=0 k = -« 


(B.24) 


3j . is defined in equation (A. 28). S. denotes the derivative with respect to its argument. 


O-iL 

\ = (A-2p) 2 + m (-^) 

. * n , n 


G 2 ‘ = (A-2p-l) 3 + nrO-^) 

n n 


(B.25) 


P # , k X! 2 ft( 2 P +1 ) ^ 1_e E X-l,p,Q-H,A+2 ' ^‘ 2p ) E A-l,p,Q+l,A+l 3 e i,o^.-(q+i),o 
1* Q= L 

t(2 " 2p+1 ^ E ^.p^-M + 2 + <«» E ^ PjQ -M +1] 

■ 6 E ^p,Q, 2+1 ,{[ l A-p+1] E l,0,k-(Q+1),0 ■ [ l Ahp+l1 E -l J 0 J k.(Q-l),0 5 

+ [E ^+ 1 ,P,Q-1,^ " E A-I,p,Q+l,A+2 ] + (A " 2p) ^PjQjW 

f 2 E +1 -® ^lyO^-CQ+l) ,0 ' B -l,0,k-te+l),0 ] ,] 1 


0,0,k-Q,-l 2^^5 [(A-2p-tQ) + G.] 


[A-2p*Q) + G. 
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( ) V — ^ v"' 1 

F f p (e) . = V" Z tE A,p,Q-l,^+l E -l, 0 ,k-(Q-l),-l ' V2,p,Q+M+l E 1,0^-(Q+1),-1 } [ U“ 2 p-l-KJj-KJ p J (B.27) 


Q= -« 


P (L) (e) = £ 

q = _ c 


\p,Q,A+l • 3 ^ 1 | 1 _p 2^ (E JS+1 , Pj Q~M+ 2 " E X-l,p > Q+l,il+2 )+(i " 2p)E X,p,Q,£+3 ]5 


E 


'0,0,k-Q,l 


+ tE l,0,k-(Q+l),-a + 2(1^2) ^ E 2,o;k-(Q+2),-l " 3 E 0,0,k-Q,-1^ [ ' (2p+1 ^ E je-l,p,Q-l,i+2 + y£^- 


U) 

I 


E A- 1 ,P,Q+M+ 1 ] 


[E 


-I,0,k-(Q-1),-1 


2 ( 1 ^) 


^ E -2,0,k-(Q-2),-l " 3 E 0,0 J k-Q^-l^^ 2 ' 4 ‘ 2p+:L ^ E £+l,p,Q-l J i+2 


U-2p) 

~2 


E £+l,p,Q-l,X+l] 


20 

* ^ X ' 2p ^ E ^p ; Q,A+l ^l,0,k-(Q+l),-l 


E -1jOj k-(Q-l),-! 



(E 2,0,k-(Q+2),-l 


E -2,0^-(Q-2),-l)} 


1 

[ (•a- 2 p-+Q)-K 5 1 ] 


(B.28) 
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APPEHDXK C 

Graphs of Inclination Function F^ 
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APPENDIX D 

Graphs of Inclination Function F 
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Graphs of Eccentricity Function 
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